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Abstract 



Properties of the three-nucleon bound state are examined in the Faddeev 
formalism, in which the quark-model nucleon-nucleon interaction is explicitly 
incorporated to calculate the off-shell T-matrix. The most recent version, fss2, 
of the Kyoto-Niigata quark-model potential yields the ground-state energy 
-E(^H) = —8.514 MeV in the 34 channel calculation, when the np interaction 
is used for the nucleon-nucleon interaction. The charge root mean square 
radii of the and ^He are 1.72 fm and 1.90 fm, respectively, including the 
finite size correction of the nucleons. These values are the closest to the 
experiments among many results obtained by detailed Faddeev calculations 
employing modern realistic nucleon-nucleon interaction models. 



21A5.+V, 21.30.-x, 13.75.Cs, 12.39.Jh 



1 



Typeset using REVT^ 



I. INTRODUCTION 



All the present-day quark-model descriptions of nucleon-nucleon {NN) and hyperon- 
nucleon (YN) interactions incorporate important roles of the quark-gluon degrees of free- 
dom in the short-range region and the meson-exchange processes dominating in the medium- 
and long-range parts of the interaction. For example, the Kyoto-Niigata quark-model po- 
tential employs a one-gluon exchange Fermi-Breit interaction and effective meson-exchange 
potentials (EMEP's) acting between quarks, and has achieved accurate descriptions of the 
A^A^ and YN interactions with a limited number of parameters. The early version, 

the model FSS includes only the scalar (S) and pseudoscalar (PS) meson-exchange 

potentials as the EMEP's. This model is superseded by the new model fss2 |Q, which has 
also introduced the vector (V) meson-exchange potentials and the momentum-dependent 
Bryan-Scott terms included in the S and V meson EMEP's. Owing to the introduction of 
the V mesons, the model fss2 in the A^A^ sector has attained the accuracy comparable to 
that of one-boson exchange potential (OBEP) models. For example, the values defined 

by = J2iLi [^i'^'' ~ S^^^^ /N for the J < 2 np phase-shift parameters in the energy range 
^lab = 25 - 300 MeV are a/x^ = 0.59° for fss2, which is compared with the values 1.10°, 
1.40° and 1.32° for OBEP, Paris and Bonn potentials, respectively. The incorporation of 
the momentum- dependent Bryan-Scott term is favorable in extending our quark-model de- 
scription of the A^A^ scattering at the non-relativistic energies to the higher energies up to 
Tiah = 800 MeV, and also in describing reasonable asymptotic behavior of the nucleon s.p. 
potentials in the high-momentum region. The agreement of the higher partial waves up to 
J=4 with the phase shift analysis is also improved. In both models FSS and fss2, the exist- 
ing data for the YN scattering are well reproduced and the essential feature of the AN-HN 
coupling remains almost unchanged. Fixing the model parameters in the strangeness S = 
and —1 sectors, we proceed to explore interactions for any arbitrary combinations of octet 
baryons [Bg). @] The BsBg interactions in 5* = —2, — 3 and —4 sectors include the AA 
and SN interactions, which are recently attracting much interest in the rapidly developing 
field of the hypernuclei and the strangeness nuclear matter. The interaction derived in these 
models may be used for realistic calculations in few-baryon systems, like the triton ^H, the 
hypertriton ^H, and also in various types of baryonic matter. This project, however, involves 
a non-trivial problem of determining how to extract the effective two-baryon interaction from 
the microscopic quark-exchange kernel. The basic baryon-baryon interaction is formulated 
as a composite-particle interaction in the framework of the resonating-group method (RGM). 
If we rewrite the RGM equation in the form of a Schrodinger-type equation, the interaction 
term becomes non-local and energy dependent. Furthermore, the RGM equation sometimes 
involves redundant components, due to the effect of the antisymmetrization, which is related 
to the existence of the Pauli-forbidden states. In such a case, the full off-shell T-matrix is not 
well defined in the standard procedure, which usually assumes simple energy-independent 
local potentials. Since these features are related to the description of the short-range part 
in the quark model, it would be desirable if the quark-exchange kernel could be used directly 
in application to many-baryon systems. 

In this paper, we will show some results of the Faddeev calculation which directly employs 
the quark-model A^A^ interactions fss2 and FSS to derive the off-shell T-matrix. Following 
the notation in ref. we write the RGM equation of the (3g)-(3g) system in the form of 



2 



the Schrodinger-type equation 



X = 0, (1.1) 



where e is the total energy in the center-of-mass system, measured from the two-cluster 
threshold, e = E — 2En, ho is the kinetic-energy operator of the A^A^ relative motion, and 

V^^^^ie) = VY, + G + eK , (1.2) 

is the RGM kernel composed of the direct potential Vd, the sum of the exchange kinetic- 
energy and interaction kernels, G = + G^ , and the exchange normalization kernel K. 
Since there is no Pauli forbidden state in the NN system, we can solve the Lippmann- 
Schwinger equation 

Tico,e) = V^^«^(e) + V^^«^(e)^?o(c^)r(u;, e) , (1.3) 

with go{uj) = l/{uj — ho + iO), by assuming the 2N energy e as a mere parameter. The 
Faddeev equation for the 3A^ bound state is given by the eigenvalue problem 

A(£;)V„ = Gi-'\E)Ti'\E,eM(3 + i^^) , (1.4) 

with X{E) = 1, where the two-body T- matrix in the three-body model space is given by 

Tj^\E,e^) = To,iE - hoa,ea) , (1.5) 

and Go{E) = 1/{E — Hq + iO) is the free Green's function for the three-body kinetic-energy 
operator Hq = hoa + hoa- The energy dependence of the two-cluster RGM kernel is self- 
consistently determined through 

= {^\hoo. + V^'^^ie^m = ^E + ^{^M^ , (1-6) 

where = (/Jq, + + yj^ is the normalized total wave function for the 3A^ bound state. 
In practice, we start from some specific values of Ea and E, and solve Eq. ( |1.4| ) to find a 
negative three-body energy E such that the eigenvalue X{E) becomes 1. The normalized 
Faddeev component ipa yields a new value of Ea through Eq. (|1.6|) .p| Since it is usually not 
equal to the starting value, we repeat the process by using the new value. This process 
of double iteration converges very fast if the starting values of Ea and E are appropriately 
chosen. 

For the numerical calculation, we discretize the continuous momentum variables p and 
q for the Jacobi coordinate vectors, using the Gauss- Legendre rii- and n2-point quadrature 
formulas, respectively, for each of the three intervals of - 1 fm~^, 1-3 fm~^ and 3-6 fm~^. 
The small contribution from the intermediate integral over p beyond po = Q fm~^ in the 
2N T-matrix calculation is also taken into account by using the Gauss-Legendre ns-point 



^As a system of identical three particles, £a, and e-y are all equal and are expressed as e in 
Tabled. 
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quadrature formula through the mappings = po + tan(7r/4)(l + The momentum region 
q = 6 fm~^ - oo is also discretized by the point formula just as in the p discretization case. 
We take ni-n2-n3=10-10-5, for which well converged results are obtained at least for 2 and 5 
channel calculations. The partial-wave decomposition of the 2N RGM kernel is carried out 
numerically using the Gauss-Legendre 20-point quadrature formula. The modified spline 



interpolation technique developed in |10| is employed for constructing the rearrangement 



matrix. For the diagonalization of the large non-symmetric matrix, the Arnordi-Lanczos 



algorithm recently developed in the ARPACK subroutine package |^ is very useful. 

Tables I and |T| list the deuteron properties and the NN effective range parameters pre- 
dicted by fss2 and FSS, respectively. All the calculations in the present paper are carried 
out in the isospin basis. For a realistic calculation of the binding energy, it is essen- 
tial to use the NN interaction that reproduces the correct D-state probability (Pd) of the 
deuteron and the effective range parameters of the ^Sq scattering. [|12| Since all the real- 



istic A^A^ interactions reproduce the A^A^ phase shifts more or less correctly, the strength 
of the central attraction is counterbalanced with that of the tensor force. Namely, if the 
interaction has a weaker tensor force, then it should have a stronger central attraction. Gen- 
erally speaking, the effect of the tensor force is reduced in the nuclear many-body systems, 
in comparison with the bare two-nucleon collision. This implies that the A^A^ interaction 
with a weaker tensor force is favorable, in order to obtain sufficient binding energies of the 
nuclear many-body systems. The weak tensor force, however, causes various problems like 
a too small value for the deuteron quadrupole moment Qd and some disagreement of the 
mixing parameter ei of the ^S*! + ^Di coupling. For example, the Reid soft core potential 
(RSC) [|I^ gives Pd = 6.5 % and predicts too small binding energy, Bt = 7.35 MeV, 
compared with the experimental value B^^^ = 8.48 MeV. A series of the Bonn potentials 
reproduce the A^A^ phase shifts very accurately, but they have a tendency that the tensor 
force is generally rather weak. The model C has the strongest tensor force Pd = 5.61 %, 
yielding Bt = 7.99 MeV. The value Pd becomes smaller for models B and A, and the value 
of Bt becomes larger, correspondingly. The following results are given in ref. 0]: model-B 
{Pd = 5.0 %, Bt = 8.13 MeV), model-A {Pd = 4.4 %, Bt = 8.32 MeV). These results 
are all obtained in the 34 channel calculations (including the 2A^ total angular momentum 
J < 4), and by using the np interaction. In fact, the effects of the charge dependence and the 
charge asymmetry are important for the detailed discussion, and it is estimated to be about 



190 keV in refs. and [|14[. The most recent Faddeev calculation employing the CD-Bonn 
potential [|r^ incorporates these effect, and predicts Bt = 8.014 MeV for Pd = 4.85 %. 
The present status of the '^H binding energy calculation is summarized as that more than 0.5 
MeV is missing if the two-nucleon force of any realistic A^A^ interactions is only employed. 



m 



On the other hand, our result of Pd in Table| is about 5.5 % both in the fss2 and FSS 
cases. We think that this is a reasonable value, in spite of the fact that Qd of fss2 is too 
small. This is because a careful evaluation of the meson-exchange current contributions to 
Qd, which could be as large as 0.01 fm^ [|T3,0], must be made. Our results of the effective 



^ These na points for p are not included for solving the Faddeev equation (0), since it causes a 
numerical inaccuracy for the interpolation. 
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range parameters in Tables|| are not as perfect as those of the Bonn B potentiaL It should 
be noted that the effects of the higher-order terms of the Coulomb interaction are not 
incorporated in these calculations. The deuteron binding energy and the scattering length 
tts for the ^5*0 state are fit in determining our quark- model parameters. 

Table fT\\ list the results of the Faddeev calculations by fss2 and FSS in various types of 
truncations of the model space. The 5 channel calculation with J < 1^ incorporates only 
the partial waves ^5*1 + ^Di and ^5*0 for the 2N T-matrix. Similarly, the 18 and 34 channel 
calculations incorporate the partial waves with J < 2 and J < 4, respectively. We find that 
the energy gain in the 5 channel to 34 channel calculation is about 330 - 360 keV, which is the 
same tendency for the realistic A^A^ potentials with a strong tensor force, such as the RSC 
and Paris potentials [l^. The convergence is not enough even in the 34 channel calculation, 



and we expect the further energy gain of the order of 10 - 20 keV. The model fss2 predicts 
Bt = 8.51 MeV and seems to give too large binding energy, compared with experiment. In 
fact, it underbinds by 150 - 160 keV, if the effects of the charge dependence and the charge 
asymmetry of the A^A^ interaction is taken into account. The scenario assuming the most 
favorable Bonn A potential is given in Table 11.1 of ref. 0, which tells us that the corrected 
value due to the charge dependence and the charge asymmetry of the two-body force is 
8.13 MeV and the rest, 350 keV, is attributed to the combined contribution of the three- 
body force and the medium effect of the two-body force. Our result using the quark-model 
potentials indicates that one can reduce the net effect besides the two-nucleon force to less 
than half of the OBEP values, by keeping the deuteron D-state probability in a reasonable 
magnitude. 

Note that the 2A^ energy Sa in Eq. ( |1.6| ) is directly related to the separation of the total 
energy E(^B.) into the kinetic-energy contribution, {Hq) = 2(35q, — E), and the potential- 
energy contribution, {V) = 3{E — 2ea)- In the 34 channel calculations, these are given by 
{Ho) = 43.95 MeV, {V) = -52.47 MeV for fss2, and (Ho) = 41.83 MeV, {V) = -50.22 MeV 
for FSS. If we compare these with the results of the CD-Bonn potential ((-f^o) = 37.42 
MeV, {V) = -45.43 MeV) and the AV-18 potential {{Ho) = 46.73 MeV, {V) = -54.35 
MeV), we find that our quark- model potentials give a moderate amount of the kinetic-energy 
contribution just between the CD-Bonn and AV18 potentials. 

Table PI also shows the calculated charge root mean square (rms) radii of and ^He 
obtained by fss2 and FSS. The finite size corrections of the nucleons are made through 

(r2)3H = i?c(^H)^ + (0.8502)2 _ 2 ^ (0.3563)^ , 

(r2)3He = RcC'^ef + (0.8502)2 _ 1 ^ (0.3563)^ , (1.7) 

where Rc^ stands for the square of the charge rms radius for the point nucleons. Since our 
3A^ bound state wave functions are given in the momentum representation, we first calculate 
the charge form factors Fc{Q'^), according to the formulation given in ref. |2^. Rc^ is then 



extracted from the power series expansion of Fc{Q'^) with respect to Q^. We have employed 
20 points, Q = 0.05 x n fm~^ with n = 1 - 20, for the extrapolation to Q = 0. In the present 
calculation, the Coulomb force and the relativistic correction terms of the charge current 



operator are entirely neglected. The experimental values are rather difficult to determine, 
as discussed in ref. [^]. Here we compare our results with two empirical values 
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1.70 ± 0.05 fm [26] rjYs — _ j 1.87 ± 0.05 fm [ge 



,^^)^H { 1 oi _i_ n riK IFvtT ' \ ('"^)3He { i no _i_ n no [Fv^l • (l-S) 



1.81 ± 0.05 f m ' V \' /^He | ^ gg ^ q 

We find that the agreement with the experiment is satisfactory both for fss2 and FSS. 

The Faddeev calculations for ^H, using the quark-model A^A^ potentials, have been carried 
out by Takeuchi, Cheon and Redish p8|, and recently by the Salamanca- Jiihch group p9 



In the former calculation, the model QCM-A, by the Tokyo University group gives the A^A^ 
phase shifts with almost the same accuracy as our model FSS. The model QCM-A predicts 
Pjj = 5.58 % for the deuteron D-state probability and Bt = 8.01 - 8.02 MeV for the 
binding energy in the 5 channel calculation. This is very similar to our results for the model 
FSS. On the other hand, the Salamanca- Jiilich group predicts Bt = 7.72 MeV, in spite of 
the very small Z)-state probability Pd = 4.85 %. It is not clear to us how they treated the 
energy dependence of the RGM kernel at the process of the separable expansion for solving 
the Faddeev equation. They have to improve the fit of the A^A^ phase shifts for higher 
partial waves (especially, for the P waves), in order to extend their calculation to more than 
5 channels. 

In summary, we have carried out the Faddeev calculation for the three-nucleon bound 
state, by explicitly incorporating the off-shell T-matrix derived from the RGM exchange 
kernel of the quark-model A^A^ interaction. The energy dependence of the two-cluster RGM 
kernel is self-consistently treated. For the two models fss2 and FSS @J^, we 

have obtained ^(^H) = -8.514 MeV (fss2), and -8.390 MeV (FSS) in the 34 channel 
calculation using the np interaction. The charge rms radii of the and ^He are in fair 
agreement with experiment: ^ {r'^)^ = 1.72 fm (fss2), 1.74 fm (FSS) and -^(7^^)311^ = 
1.90 fm (fss2), 1.92 fm (FSS). In these calculations, the Coulomb force and the relativistic 
correction terms are neglected. In view of the fact that the A^A^ phase shifts of FSS are not 
that excellent, the results of fss2 are more meaningful. These numbers are the closest to 
the experiments among many results obtained by Faddeev calculations employing modern 
realistic A^A^ interaction models. Since both models fss2 and FSS have a common feature 
in describing the short range correlation by the quark exchange kernel, it is important to 
clarify the mechanism in which the quark-model potentials give larger binding energy 
than the meson-exchange potentials. The off-shell behavior of the RGM T-matrix is closely 
connected to this alternative description of the short-range correlations. More detailed study 
on this point is now under way. 
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TABLES 



TABLE L The deuteron properties by fss2 and FSS in the isospin basis. The results by the 
Bonn B potential |^ are also shown for comparison. A small difference in FSS from Table IV of 
is due to the numerical inaccuracy in the previous calculation. The effect of the meson exchange 
current is not included in the calculated values of Qd and fid- 





FSS 


fss2 


Bonn B 


Expt. 


Ref. 


ed (MeV) 


2.256 


2.225 


2.2246 


2.224644 ± 0.000046 




2C| 


Pd (%) 


5.86 


5.49 


4.99 








rj = Ad /As 


0.0267 


0.0253 


0.0264 


0.0256 ± 0.0004 




211 


rms (fm) 


1.963 


1.960 


1.968 


1.9635 ±0.0046 




2C| 


Qd (fm^) 


0.283 


0.270 


0.278 


0.2860 ±0.0015 




22] 




0.8464 


0.8485 


0.8514 


0.857406 ± 0.000001 




23 1 



TABLE IL The A^A^ effective range parameters calculated by fss2 and FSS in the isospin basis. 
The results by the Bonn B potential [0 are also shown for comparison. The higher-order terms of 



the Coulomb force are not included. The experimental values are taken from [20|. 



FSS fss2 Bonn B Expt. 

as (fm) -23.64 -23.76 -23.75 -23.748 ± 0.010 

(fm) 2.62 2.58 2.71 2.75 ±0.05 

as (fm) 5.41 5.399 5.424 5.424 ± 0.004 

rs (fm) 1.76 1.730 1.761 1.759 ± 0.005 
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TABLE III. The three-nucleon bound state properties predicted by the Faddeev calculation 
with fss2 and FSS. The np interaction is used in the isospin basis. The discretization points of p and 
q are specified by the values of ni-n2-n3=10-10-5 (see the text). The column "channels" implies 
the number of two-nucleon channels included, and n-max = n(3ni)(3n2 + ns) for the n channel 
calculation is the dimension of the diagonalization for the Faddeev equation. E{^}1) is the ground 
state energy, and VTr^JsH and V ('^^)3He are the charge rms radii for and ^He, respectively, 
with the proton and neutron size corrections introduced by Eq. (1/7). The Coulomb force and 
the relativistic corrections are neglected, e is the 2N expectation value, Eq. ([L.6|), determined 
self-consistently. 



model 


channels 


^max 


ei2N) 














(MeV) 


(MeV) 


(fm) 


(fm) 




2 ch 


2,100 


2.361 


-7.807 


1.80 


1.96 




5 ch 


5,250 


4.341 


-8.189 


1.75 


1.92 


fss2 


10 ch 


10,500 


4.249 


-8.017 


1.76 


1.94 




18 ch 


18,900 


4.460 


-8.439 


1.72 


1.90 




34 ch 


35,700 


4.488 


-8.514 


1.72 


1.90 




2 ch 


2,100 


2.038 


-7.674 


1.83 


1.99 




5 ch 


5,250 


3.999 


-8.034 


1.78 


1.95 


FSS 


10 ch 


10,500 


3.934 


-7.909 


1.78 


1.97 




18 ch 


18,900 


4.160 


-8.342 


1.74 


1.93 




34 ch 


35,700 


4.175 


-8.390 


1.74 


1.92 
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